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arising in the study of symplectic ensembles of random matrices. 

PACS numbers: 02.30.Gp, 05.45.Mt 



! 1. INTRODUCTION 

o 



o 



After the publication of (l| and 0, late Prof. M. L. Mehta was not particularly happy with the formal nature of 
the presentation. Following his suggestion, this paper attempts to give an explicit workout of the skew-orthogonal 
polynomials (SOP) corresponding to the quartic potential. 

We study SOP arising in the study of symplectic ensembles of random matrices. The corresponding weight function 
is w{x) = exp[— y(a;)], where 



vix) = - + —, aeM, (1.1) 

J_i ' is the quartic potential. In this paper, we outline an iterative technique to develop these polynomials which can be 
used to obtain the level-density and 2-point function for the symplectic ensembles of random matrices. Here, we 



■ emphasize that this method can be easily extended to all potentials of the form 

n-) = E^^, «2. = i. (1.2) 

k=l 

> 

C — ■ Without loss of generality, we will be dealing with monic SOP of the form 

\o ' 

O; (/.„(x)=w;(a;)E4"^x^ 4") = 1, (1.3) 

o 

■ They satisfy the skew-orthonormalization relation 



and define 



(j)n{x)^m{x)dx ^ gnZnrn, (1-5) 



Z = [ ^,1 )+...+, (1.6) 



where 

. - f 

is an anti-symmetric block-diagonal matrix with = —1. (7„ is the normalization constant with the property 

g2n^ 9271+1- (1-7) 

Here ipn(x) is a polynomial of order n + 2d ^ 1 with leading coefficient —1. We have dropped the superscript /3 (as in 
Q) and H) since we will only be interested in SOP arising in the study of symplectic ensembles of random matrices. 
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As explained in Ref. [l[ and 0], we expand ip{x) in terms of (f>{x) and write 

A; 

X1pn{x) = ^ Rn,n<Pni{x), (1.8) 
m=j 

where from Eqs. ()1.2l 11.31 ll-4|) . one can say that k — n + 2d. Furthermore, using the integral ^^xipn{x)ipm{x)dx, we 
can show that the matrix R satisfy anti-self dual relation 

R = ZR*Z=-R°. (1.9) 
Summing up these results, we can write for any polynomial weight 

Xil2n{x) = i?2n,2n+2d</>2n+2d(a;) + . . . + i?2ri,2n-2d02ri-2d(a;), (1-10) 
X4'2n+l{x) = R2n+l,2n+2d+l4'2n+2d+lix) + ■ ■ ■ + R2n+l,2n-2d4'2n-2dix), (l-H) 

where (/>_m = 0, m being a positive integer. Our choice of the potential {u2d = 1) ensures that 

Rm.,m+2d = -l, m>0. (1.12) 

Also from Eq. (|1.9p . we can show that 



1, 


n > d. 


0, 


n < d, 


1, 


n > d, 


0, 


n < d. 



(1.13) 

Following this brief recapitulation of [l| and 0], we will outline our plan of action. 

1. We will calculate the SOP's (j>n{x) for < n < 2d (for the particular case d = 2) using generalized Gram-Schmidt 
orthogonalization. 

2. Using these polynomials, we will use Eqs. (|1.10|) and (jl.lip to calculate the higher order polynomials recursively. 
The coefficients will be expressed in terms of integrals of the form x" exp[—2V {x)]dx . 

3. We show that these coefficients themselves satisfy a set of difference equations. 

4. Using these difference equations, we also obtain the normalization constants. 

5. Finally, we talk briefly about the zeros of these polynomials. 



2. THE GRAM-SCHMIDT TECHNIQUE 

Putting 71 = in Eqs. (|1.10p and p. lip (keeping in mind Eq. (jl.l3p ). we can see that it is useless, unless we have 
information about the SOP's (j)m for 2d > m > 0. To overcome this problem, we will obtain these 2d polynomials 
{(j)o{x), . . . , (j)2d-iix)) using the Gram-Schmidt technique for SOP. From here onward, we will focus our attention on 
the specific weight function defined in Eq. p.ip . although the method outlined can be extended to any d. For d = 2, 
the first 4 monic polynomials can be written as 

(poix) = w{x), (j)i{x) ^ xw{x), (2.1) 

02(x) = {x^ + c''o^)w{x), Mx) = (x^ + c^^^x)w(x). (2.2) 

Correspondingly 

^x) = -V'{x)w{x), iPi{x) = {I - xV'{x))w{x), (2.3) 

^2{x) = [2x ^ {x^ +c'a^)V'{x)]w{x), ipsix) = [Sx"^ + cf^ - V ix){x^ + cf ^ x)]w{x). (2.4) 

Using 



/oo 
(j)2{x)^pi{x)dx — 0, 
-oo 



(2.5) 
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we get 



Similarly, using 



we get 



We also get 



° J^^{l-xV'{x))w^x)dx ■ ^ ' ' 



Mx)Mx)dx = 0, (2.7) 



[°° x^V'(x)w'^(x)dx 
^ J^^xV'{x)w^{x)dx 



oo 



50 = 51=/ (f>o{x)ilJi{x)dx, 52=ff3=/ (?!)2(a;)'03(a;)rfa;- (2.9) 

J — OO 'J — oo 

3. RECURSION RELATION 

Now, let us look at Eqs. lfTTO)) and (frTTj) . Since V{x) = V{-x), we have 

''p2nix) -'02n(-a;), ■02n+l(2^) = 1p2n+li-x). (3.1) 

Thus the odd (even) terms will be absent in Eqs. p.lO|) ( (|1.11|) ) since 



oo 



-R2n,2n+2*:+l = / xij2n{x)lp2n+2k{x)dx = 0, 

52n+2fc J-oo 

and 

1 /"^ 

^2n+1.2n+2A; = / X'llj2n+l{x)lp2n+2k+l{x)dx = 0. (3.2) 

g2n+2k J-oo 

Using these results and Eq. (|1.12p . we can rewrite Eqs. (|1.10|1 and (jl.lip for n = and d = 2 as 

x^l^aix) = -(j)4{x) + Ro^2(l)2{x) + Rofl(f)o{x), (3.3) 

x^piix) = -05(a:) +i?i^3(/)3(x) + i?i,i0i(a;). (3.4) 

Also, using the skew-orthogonality property, we get 

1 f°° 

Rofl = — / xipo{x)'ilJi{x)dx, (3.5) 



50 J-c 
1 

52 



i?o,2 = — / a:V'o(a;)^3(a;)rfa;- (3.6) 



Similarly, 



1 f°° 

Ri,i = / x'il:oix)'ipi{x)dx ^ -Rq^o, (3.7) 

90 J-oo 

1 r°° 

Ri.3 = / x'ilJi{x)tl;2ix)dx. (3.8) 

52 J-oo 

Here, we note that Rq q — — can be obtained directly from Ea. (|1.9p . Once the coefficients are known, one can 
calculate ^4(0;) and 05 (a;) using Eqs. (|3.3|) and (|3.4|) . Using p.4p . we can calculate tp4{x) and i^six). We can also 
calculate 

/•oo 

54 = 55= / (^4(a;)V'5(a;)da;. (3.9) 
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54 J-oo 
1 

94 



R3.5 = — - I xtp3{x)ilj4{x)dx, (3.14) 



00 



We will now calculate (j)e{x) and (j)7{x) using 

xip2{x) = -06(2;) + i?2,4(/'4(a:) + i?2,202(a;) + i?2.o?!'o(a;), (3.10) 

Xipsix) = -(j>j{x) + R3,5(j>i,{x) + R3,3<j)3{x) + R3^i(l)i{x). (3.11) 

Again, using Eq. ()1.9p . we have 

-^2,0 — Rs.l — ~-Ro,2, -^2,2 = — -^3,3- (3-12) 

i?o,2 and i?i,3 has already been calculated in (|3.6p and (|3.8p respectively. Also 

1 /""^ 

i?2,4 = — / X'ip2{x)ij5{x)dx, (3.13) 



and 

1 f"" 

i?2.2 = -i?3,3 = / X'lp2{x)i^3{x)dx. (3.15) 

92 J-oo 

With these, we can calculate (/'6(a^), 4'i{x) and correspondingly "06 (2:) and ip7{x) and have 

/oo 
(/)6(a;)V'7(a;)da;. (3.16) 
-C30 

Following the same technique, we can obtain the polynomials for all n. For n > 2 and d = 2, we may rewrite Eas. (|1.10p 
and pTTTj) as 

2 2 

X^27i{x) = ^ i?2n,2n+2m<?!>2n+2m(a;), a;'(/'2Ti + l (a^) = ^ -R2ri+l,2n+2m+l 02ri+2m+l (a:) . (3-17) 

m— — 2 m— — 2 

From Ea. (|1.12p . we have 

^2n,2n+4 ~ -R2n+l,2n+5 = ^1- (3.18) 

Also from Eq. (|1.13|) . we have 

-R2n,2n-4 = ^-R2n-3,2n+l = Ij ^2ri+l,2n-3 = ^^2ri-4,2n = 1- (3.19) 

R2n,2n-2 (for the evcn case) and i?2n+i,2ri-i (for the odd case) is already known since 

R2n,2n-2 — ^ ^^2n-l,2ri+l , -R2n+l,2ri-l = ~^2n-2,2n- (3.20) 

Finally, we are left with terms of the form Rk.k+2 and Rk,k (for both k odd and even). They are given by 

1 f°° 

R2n,2n+2k = / X1p2n{x)'lp2n+2k+l{x)dx, fc = 0, 1 (3.21) 

92k+2n J-oo 



oo 



-R2n+l,2n+2fc+l = / X^2n+l{x)llj2n+2k{x)dx , k = 0,l. (3.22) 

g2k+2n J-oo 

Here, we note that by expanding xij}2n (or a;V'2n+i), we can evaluate i?2n,2n+2fe (or i?2ri+i,2n+2/c+i) analytically. But 

that will involve terms like <^2n+2'' (^r '^2ri+3^^ ) thereby cannot be used to obtain the polynomials recursively. We 
might recall that for monic orthogonal polynomials, this problem does not exist However it is possible to evaluate 
Eqs. (|3.2ip and (|3.22p numerically, although it may be a tiresome process. 
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4. THE NORMALIZATION CONSTANT 

So far, we have outlined a formalism to obtain the polynomials recursively, where the recursion coefficients are 
expressed in terms of certain integrals which needs to be evaluated at every iteration. It is not practical to use this 
process for the study of large n behavior of these polynomials. Nor is it convenient to study potentials with larger d, 
since the number of terms in the recursion relation increases with d. 

In this section, we present an alternative technique to obtain both the recursion coefficients and the normalization 
constant. We begin with the identity (obtained by integration by parts) 

j [x{x%lj-j{x))"\ (l)k{x)dx = - J {x^pkix)){x^pj{x))'dx, j, /c = 0, 1, . . . . (4.1) 

The only criteria to use this formalism is that as initial condition, we need to know these polynomials for n — 
0, . . . , 2(i — 1. Taking j — 2n and fc = 2n + 1 in Eq. (|4.ip . we get the recursion relation for the normalization constant. 
We get 

54+7o52-(l + 7o)go = 0, n = Q, (4.2) 
56 + 7iff4 - (1 + 70 +7i).92 +7030 = 0, n = 1, (4.3) 

ff2n+4 + 7n52n+2 - (2 + 7n + 7n-l)52n + 7n-lff2n-2 + ff2n-4 = 0, 71 > 2, (4.4) 

where 

7„ := i?2n,2n+2^2n+l,2n+3- (4-5) 

Here, we have expanded xip2n{x) and xip2n+i{x) (using Ea. (|3.17p ). and using Eqs. (|3.18p . (|3.19p and (|3.20p '). get the 
result. We have assumed (p-n and consequently g-n and 7_„ is zero, n being a positive integer. Here, we must also 
remember that for these SOP, g2n = 32n+i- 

To evaluate 7, we need to know Rk,k+2 for both k odd and even. This can be calculated from Eq. (|4.ip by putting 
j = 2n, k = 2n + 3 and j = 2n + 1, fc = 2n + 2. We get 

R2n+3,2n+5{g2n+4 — 32ri+2) ^ ^2n+l,2n+1^2n,2n+2 (ff2n — 92n+2) + ^2ri- l,2ri+l (32n-2 — 52n+2) = 0, (4.6) 
-R2n+2,2n+4(ff2n+4 ^ 32ri+2) — ^2n,2n^2n+l,2n+3(ff2n " 32n+2) + ^2n-2,2ri(32n-2 " ff2n+2) = 0. (4.7) 

These three equations can be used together to obtain gk for all k. For example, a knowledge of i?o.2 and -R1.3 (which 
we get from Eqs. (l3.6p and p.Sp ) will give 79. This will give 54 from (|4.2p (again we need g^ and 32 which we will 
extract from Eq. (|2.9p V A knowledge of 34 in turn yields i?3^5 and i?2,4 (where we use _Ro.o a-nd -Ri,i, calculated from 
Eq. p.Tp ) from Eqs. (|4.6p and (|4.7p . which gives 71. This in turn can be used to calculate g^ ()4.3p and so on. 

However, we still need to know Rj^j for j > 2. This can be obtained by writing j = 2ri + 3 and fc = 2ri in Eq. (|4.ip . 
We get 

^2n+3,2n+3^2 n,2n+2 \y2n-\-2 - g2n) + R2 n+3,2n+5 \92n ^2n+4 ) + R 2n-l,2n+l(ff2n " 52n-2) — 0. (4.8) 

^2n+2,2n+2 Can be calculated from the relation i?2n+2,2n+2 = — ^2n+3,2n+3- Thus a knowledge of (74, i?o,o and i?o,2 
gives i?3^5 (|4.6p . Then using 34, 52, ^3,5 and i?o,2, we can obtain i?3_3 (and hence -R2,2) from (|4.8p . 

To summarize, we have outlined a recursive technique to obtain the normalization constant. This necessitates a 
knowledge of coefficients of the form Rk,k+2- Again to evaluate this (also needed in Eq. ()3.17p ). we need Rk,k- We 
obtain a series of self-consistent recursion relations to obtain these coefficients. 

Comment: In the context of random matrix theory, we may point out that now we can obtain various statistical 
properties like the level-density and 2-point functions for even finite dimensional symplectic ensembles of random 
matrices (jH (3). 

5. THE ZEROS OF THESE POLYNOMIALS 

Having obtained the polynomials for the specific case of d = 2, we will discuss briefly about the zeros of these SOP 
in this section. We know that in a given interval [a, b] 



6 



VriT^O, / ipn{x)w{x)dx = 0, n>2. (5.1) 

J a 

So ^pn{x) should have atleast one point in the interior of [a,b] where it changes sign. Let there be k such points 
Xi,X2, ■ ■ ■ ,Xk. Then the function ijjn{x){x — xi) . . . (x — Xf^) is positive or negetive definite for k < n + 2d — 1, since 
ipnix) is a polynomial of order n + 2d — 1. This implies 

b 

ipn{x){x — xi) . . . {x — Xk)dx ^ 0, yk < n + 2d ~ 1. (5-2) 

However, from skew-normalization relation, this condition is satisfied if and only if 

1. For n — 2m, k = 2m + 1. 

2. For n = 2m +1, k — 2m. 

3. {x - xi) . . .{x - Xk) = (f)kix). 
This implies that 

a. (pkix) is a polynomial of order k with k real zeros. 

b. ipk{x) is a polynomial of order fc + 2(i — 1, with ^2m{x) having 2m + 1 and ^2m+i{x) having 2m real zeros. 



6. CONCLUSION 



This paper gives a formalism to derive the SOP arising in the study of symplectic ensembles of random matrices 
. Here, we would like to point out that this paper contains almost no new results which were not present in [l| and 
. However, a lot of the properties were overlooked. For example, the explicit form of the recursion coefficients, the 
necessity to use Gram-Schmidt method to obtain polynomials of order n < 2d, the zeros of 4in{x) and ipn{x) etc. At 
this point, one can easily obtain these polynomials, although a deeper insight into the recursion coefficients (specially 
its large n behavior) is an absolute necessity. However, with all its limitations, we hope that this paper will take us 
further towards our goal, which is "to develop the theory of skew orthogonal polynomials until it becomes a working 
tool as handy as the existing theory of orthogonal polynomials" ?] . 
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